Thermometry and signatures of strong correlations from Raman spectroscopy of 

fermionic atoms in optical lattices 
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We propose a method to directly measure the temperature of a gas of weakly interacting fermionic 
atoms loaded into an optical lattice. This technique relies on Raman spectroscopy and is applicable 
to experimentally relevant temperature regimes. Additionally, we show that a similar spectroscopy 
scheme can be used to obtain information on the quasiparticle properties and Hubbard bands of 
the metallic and Mott-insulating states of interacting fermionic spin mixtures. These two meth- 
ods provide experimentalists with novel probes to accurately characterize fermionic quantum gases 
confined to optical lattices. 
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I. INTRODUCTION 



Fermionic ultracold atom physics has witnessed un- 
precedented experimental progress since a quantum de- 
generate Fermi gas was first prepared in a three di- 
mensional optical lattice The recent evidence for a 
fermionic Mott insulating state @, H[ serves as a clear ex- 
ample of these rapid advances. However, despite all these 
breakthroughs, conducting experiments with fermionic 
atoms still remains a major challenge as very few probes 
are available to accurately characterize these systems. 
Among all difficulties encountered by experimentalists, 
the lack of reliable methods to adequately measure the 
temperature of Fermi gases confined to optical lattices is 
often cited as one of the main obstacles. 

Although this problem is quite persistent when con- 
fronted to lattice systems, the situation is much better in 
the continuum. There, several techniques were success- 
fully implemented to estimate the temperature of both 
fermionic and bosonic quantum gases [1, H, @|. In the 
presence of an optical lattice, different schemes to de- 
termine the temperature have been devised and experi- 
mentally tested for bosonic atoms. One approach relies 
on the direct comparison of experimental and theoreti- 
cal time-of-ffight images obtained from computationally 
expensive simulations 0- In a second method, the tem- 
perature is estimated from a measurement of the width 
of the transition layer between two spin domains cre- 
ated by the application of a magnetic field gradient Q. 
Relying on the good local resolution attainable in two 
dimensional systems, a third approach extracts temper- 
ature from the density and density fluctuations [fjj . This 
method, based on a generalized version of the fluctuation- 
dissipation theorem was analyzed theoretically in [Tol ] . 
Finally, detecting the temperature using bosonic impu- 
rity atoms which are insensitive to the optical lattice po- 
tential was put forward in pd| . 



The situation is much more difficult for fermionic 
atoms loaded into an optical lattice. In current experi- 
ments, temperature measurements are usually performed 
before switching on the lattice potential and after switch- 
ing it off [2|, |3| . However, since the temperature changes 
during the loading process, detecting the system tem- 
perature with the optical lattice on is of the utmost im- 
portance. Experimental attempts [l2| at evaluating the 
temperature of fcrmions inside a lattice were based on the 
determination of the number of doubly occupied sites. In 
this scheme, an accurate evaluation of the temperature 
requires a full theoretical understanding of the strong 
dependence of the number of doubly occupied sites with 
temperature [H, [13, [l]|. In addition, a precise exper- 
imental knowledge of the interaction strength, hopping 
amplitude, trapping configuration and particle number 
is needed. The combination of all these requirements 
renders this approach difficult to use. On the theoret- 
ical side, other methods were proposed. For example, 
one could envisage, as proposed in [nj, to extract tem- 
perature from a generalized version of the fluctutation- 
dissipation theorem used in conjunction with the knowl- 
edge of both spatially resolved system density and density 
fluctuations. However, this approach needs sufficiently 
strong density fluctuations as well as very good local res- 
olution, an experimental requirement that is far from be- 
ing met for fermionic systems. Finally, in one and two 
dimensions, measuring the intensity of the light scattered 
off the atomic lattice array was proposed to detect the 
system temperature [16j . 

In this work, we propose a novel method to measure 
the temperature of fermions loaded into two or three di- 
mensional optical lattices. The approach we put forward 
relies on transferring a portion of the atoms stored into 
the optical lattice potential to a third hyperfine state us- 
ing a stimulated Raman process. Hence, the temperature 
measurement can either be done locally or globally, and 
can be used in parallel with other probes. Depending 
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on the experimental resolution, this thermometer works 
for both free and weakly interacting fcrmions, and only 
requires the knowledge of the hopping amplitude of the 
system under study. In addition, our approach can be im- 
plemented using present fermionic ultracold atom tech- 
nology. 

Measuring the temperature of fermionic gases loaded 
into optical lattices is not the only difficulty faced by ex- 
perimentalists working with cold atoms. Indeed, identify- 
ing the different strongly correlated phases that can be re- 
alized in these systems is also a demanding task. In com- 
parison to condensed matter systems, few probing tech- 
niques are available to study cold atomic systems. Spec- 
troscopic methods [H, [a], including momentum-resolved 
radio-frequency spectroscopy [6|, have been shown pre- 
viously to be very efficient in probing characteristics of 
quantum gases not subjected to optical lattice potentials. 
Recently, it has been suggested that Raman spectroscopy 
can also be used to probe the excitation spectrum of 
strongly correlated phases of Bose gases confined to opti- 
cal lattices [l?], EH and to investigate single-particle exci- 
tations in normal and superfluid phases of fermionic gases 
In this article, we show that spectroscopy can 
be used to identify various signatures of strongly corre- 
lated fermionic phases in optical lattices. Experimentally 
detectable features include the presence of quasiparticle 
peaks in weakly and strongly correlated liquids as well as 
Hubbard bands in strongly correlated liquids and Mott 
insulators. 

The rest of the article is organized as follows: in Sec. |TT] 
we define the general setup for Raman spectroscopy. In 
Sec. IIIII we present the temperature detection scheme for 
non-interacting (Sec. IIII Ap and weakly interacting (Sec. 
IIII Bp Fermi gases confined to two- and three-dimensional 
optical lattice. Finally, in Sec. IIVI we demonstrate that 
spectroscopy can also be used to identify various signa- 
tures of strongly correlated systems. 



II. SETUP AND THEORETICAL 
DESCRIPTION 

The proposed detection schemes rely on Raman spec- 
troscopy [3, H3, This probing technique consists 
in exciting with a given energy and momentum a many- 
body state formed of a mixture of two hyperfine states 
by transferring atoms to a third state. To set the ideas 
straight, we sketch the Raman process in Fig. [TJ There 
we see that atoms from hyperfine state 1 1) are transferred 
to a different hyperfine state |3) using two Raman laser 
beams with frequencies W12 and LU23 and Rabi frequencies 
Q12 and 1^23 , respectively. The frequencies U12 and UJ23 
are both detuned from their corresponding resonances to 
state 1 2) to keep this state unoccupied. During the transi- 
tion, momentum q = ki — k2 is transferred to the atoms. 
This value can be chosen within certain bounds by ad- 
justing appropriately the angles between the two Raman 
beams and the lattice axis (Fig. [5]). For example, trans- 



ferred momentum q = could be realized using coprop- 
agating Raman laser beams. While, for many current 
lattice setups, q « (-,-) or (-,-,-) could be reached 
by aligning two counterpropagating Raman lasers along 
the diagonal of the optical lattice axes. 

Experimentally, the Raman signal is measured by 
counting the number of atoms transferred to state |3). 
This signal can in principle be resolved both in fre- 
quency and momentum. For many applications, such as 
thermometry, momentum resolved measurements are not 
needed. Nevertheless, as we will show in Sec. IIVI mo- 
mentum resolution can also provide valuable additional 
information, but achieving good momentum resolution in 
optical lattice setups is experimentally demanding. 

When only a small fraction of the atoms in hyperfine 
state |1) are transferred into state |3), the Raman sig- 
nal can usually be approximated using a linear response 
expression |l9l [20j . Within local density approximation, 
the Raman transition rate is given by 

x A(k-q,£3 k -/i -^;/i r ). (1) 

In this expression, £3 k = £3,k + V3(r) and [i r = Mo - ^i( r ) 
where £3^ is the dispersion relation for the |3) state, 
Vi,3(r) are the trapping potentials felt by states |1) and 
1 3), respectively, and /i is the chemical potential in the 
center of the trap. The momentum dependent coefficient 
W£ is due to the Wannier envelope and is given by 

= I J dr w*(r) ^, k (r) e^ * | 2 (2) 
= I J dr <k(r) u 3,k- q (r) | 2 

where u>i(r) is the Wannier function for the atoms in 
state |1) while ^3 : k(r) is the Bloch function for the 
atoms in state |3) (ui,k,U3 jk are the corresponding pe- 
riodic parts of the Bloch function). np(x) = 1/(1 + 
exp(x/(fesT))) is the Fermi function. A(k, hv;[i r ) is the 
one-particle spectral function for the (|1), |1')) mixture in 
a confining potential. The local density approximation 
has been used, so that /i r is the local chemical poten- 
tial at point r. In an homogeneous system, the spectral 
function is defined as 

„ p -E,/k B T , -E f /k B T 

A(k, hv) = £ B - Y 

x |(0/|ci k |^)| 2 8{hv + E } -Ei) (3) 

where cik detroys an atom in state |1) with momen- 
tum k, Z = J2i ex P( — Ei/ksT) is the Grand-Canonical 
partition function, the sums over i and / refer to all 
the many-body states of the system and the energy 
Ei — Ei — iiNi is rescaled by the number of particles. 
Finally, fku — ft(wi2 — ^23) is the transferred energy, and 
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the Rabi frequency is f2 e (r) = f2i2(r)f22 3 (r)/A, where A 
is the detuning. Local resolution of the Raman transfer 
could be obtained by using special configuration of laser 
beams (20| . 

From Eq. [T] we see that for a given position and mo- 
mentum the Raman spectrum is obtained from the multi- 
plication of two functions. The first function is the Fermi 
factor rip which depends strongly on temperature, but is 
independent of other parameters apart from \i Q and T 
itself. As we will show in Sec. IIII1 our temperature de- 
tection scheme relies primarily on this observation. The 
second function entering Eq. [1] is the spectral function 
which depends sensitively on the state of the system. We 
will see that this limits somewhat the possibility of a uni- 
versal temperature determination, but as we explain in 
Sec. IIVI this function provides valuable information on 
the phase of the system. 

Through out the rest of this article, we describe quan- 
tum gases confined to optical lattice potentials using the 
fermionic Hubbard model [HI, [24[ : 

H = -J ^2 ( c L c r'a + h - c ) + U^n r] h ri 

(r,r')a r 

- ^ Mr n ra , (4) 

ra 

where c\ a and c ra are the creation and annihilation op- 
erators of the fermions with a — {|1), |1')}> J is the hop- 
ping matrix element, U is the on-site repulsion, fj, r is the 
local chemical potential and n ra — c\. a c ra is the num- 
ber operator on site r. (r, r') denotes neighboring lattice 
sites. 
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FIG. 1: Atomic levels involved in the Raman process. Atoms 
in hyperfine state |1) are transferred to state |3) using two Ra- 
man laser beams with frequencies W12 and W23, respectively. 
These frequencies are both detuned by A from their corre- 
sponding resonances to state |2). States |1) and |3) are sepa- 
rated in energy by e%. 



III. TEMPERATURE DETERMINATION 

In this section, we present a novel method to evaluate 
the temperature of weakly interacting fermions confined 
to optical lattices by measuring their Raman spectrum. 
We first present this method considering non-interacting 
fermions in two- and three-dimensional optical lattices. 
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FIG. 2: Geometrical configuration of the two Raman laser 
beams carrying momentum ki and k2 respectively. 

Afterwards, we demonstrate how this procedure is also 
applicable to weakly interacting fermionic gases loaded 
into optical lattices. 

A. Temperature extraction for non-interacting 
fermions 

For the case of non-interacting fermions, i.e. U = 
0, the spectral function for trapped |1) atoms is 
A(k, hv\ (Lt r ) = 5{hv + fi r — £i,k)- In our study, we will as- 
sume that both |1) and |3) atoms are trapped by the same 
harmonic potential, i.e. Vi = V3 = Vt, and are loaded 
into an optical lattice which is felt equally by both hy- 
perfine states, i.e. induces the same hopping coefficients, 
such that £i,k = £3^ — £3 = £k- Here £3 is the energy 
offset of state |3) with respect to state |1) (cf. Fig. [TJ. 
These assumptions are valid as long as we use adequate 
hyperfine states and confine the atoms into far-detuned 
optical lattices. Under these conditions, the Raman tran- 
sition rate is given by 

r J 

x n F (e^ - hjj + £ k + Vt(t) - fJ, ) 

X (5(£3-ficJ + £k-£k-q)- (5) 

This expression depends on temperature only through 
the Fermi function rip. Therefore, detecting the sys- 
tem temperature can be done reliably by fitting i? q (w) 
with a minimum of parameters. From Eq. one also 
sees that the frequency spread of the Raman signal is 
strongly dependent on the chosen transferred momen- 
tum, q. As shown on Fig. [3l for a homogeneous system, 
atoms can only be transferred from states |1) to |3) with 
energy fvjj = e% if q = 0. Hence, i? q=0 (w) is strongly 
peaked at £ 3 /fi and zero everywhere else [25[. This fea- 
ture makes it impossible to detect the system tempera- 
ture using q = as the Raman signal is too narrow in lu. 
What we need is a Raman signal that is non-zero for a 
wide range of frequencies. This is achieved if q = (- , ^ ) 
for square lattices and q = f , f ) for cubic lattices. 
As we can see on Fig. |31 for these two configurations, 
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since £k- q = — £k, atoms can in principle be transferred 
with energies ranging from £3 — 2D to £3 + 2D where 
D is the half-bandwidth. Therefore, temperature mea- 
surements can be done with this choice of transferred 
momentum. We use the notation q = 7r/a to specify 
that we use these configurations. We also note that since 
we are interested in the momentum integrated Raman 
rate, W^ _7r ^ a can be neglected without loss of generality 
for sufficiently deep lattices as this factor only affects the 
signal overall amplitude (2r| . 



As we will show this approximate expression which has 
only temperature and a prefactor as fitting parameters is 
only valid in certain parameter regimes. However, it has 
the advantage of being very simple. 

In the following, we show that the values of tempera- 
ture and central chemical potential obtained through this 
fitting procedure agree very well with the initial system 
parameters. This means that each Raman spectrum is to 
a good extend uniquely defined by its temperature and 
chemical potential both in two and three dimensions. 



q = 



q = n/a 




-3,k-q 



FIG. 3: Left: Raman transfer with q = 0, atoms can only be 
transferred from states |1) to |3) with frequency ui = £3/?!. 
Right: Raman transfer with q = n/a, atoms can in prin- 
ciple be transferred from states |1) to |3) with frequencies 
ranging from (£3 - 2D) /h to (£3 + 2D) /h where D is the half- 
bandwidth. £i,k and £3,k are the dispersions of states |1) and 
1 3), respectively, p is the chemical potential for the mixture of 
states |1) and |l'), and £3 is the energy shift of |3) compared 
to |1). 

In the remaining of this section, we demonstrate that 
Raman spectroscopy can adequately be used as a ther- 
mometer. We conduct this demonstration in two steps. 
First, using Eq. [5j we numerically simulate a Raman ex- 
periment in which the probing lasers are shone on the 
whole system of atoms confined to a lattice and trapped 
in an harmonic potential. These simulations show that 
the resulting spectra vary significantly with temperature. 
Then, as a second step, we fit the obtained signal us- 
ing two different fitting functions. The first is the full 
continuum-space expression for the Raman rate (valid 
for large systems): 



Rv/a{u) C ( HU) 

— «r~ K 7 9v U" 



dp. 



- /,M d - 2 )/ 2 



I + e -(hCu/2+ tJ .)/(k B T) 



.(6) 



In this expression, Hlu = hco — £3, g v (e) = ^ ^ k S(e — £k) 
is the density of states of the band, d is the dimension of 
the system, V its volume and C = 2ir\fl e \ 2 /h. As we use 
the local density approximation, the Raman rate per par- 
ticle only depends on the particle number and the trap- 
ping potential through the characteristic particle number 
p = N(Vt I 'D) d / 2 where N is the total number of atoms 
in the hyperfine mixture (|1), |1')). We are left with only 
three fitting parameters: the temperature, T, the chemi- 
cal potential at the center of the trap, p D , and an overall 
multiplicative factor. The second fitting function is the 
simplified expression: 



1. Non-interacting fermions in two dimensions 

Let us first look at Raman spectra for a two- 
dimensional system with fixed characteristic particle 
number p = 2 and varying temperatures. These spec- 
tra are shown in Fig. 3J From this figure, one can see 
that all spectra can be broken down into two parts. For 
large uj, each spectrum is characterized by a signal of 
large amplitude whose shape depends on the system pa- 
rameters while, for small u>, each spectrum presents a tail 
whose shape is mainly set by temperature. The peak or 
discontinuity at w = is due to the Van Hove singularity 
of the square lattice density of states. The presence of 
sharp edges at Mj/AJ = {—2, 2} is due to the abrupt ends 
of the square lattice density of states. At low tempera- 
tures, most of the spectral weight is located at <1> > 
with a sharp step around u) ps (cf. ksT/AJ = 0.1). 
As the system temperature is increased, some weight is 
transferred into the tail and the step broadens consid- 
erably. This broadening stems from the smoothening of 
the Fermi function with increasing temperature. For high 
temperatures, the left end of the spectrum becomes sharp 
as it is cut by the edge of the density of states. 
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FIG. 4: Raman spectra for three different temperatures 
ksT/AJ — 0.1,0.5,1 at the characteristic particle number 
p = 2 in a two-dimensional system. (fcsT/4J) s is the ex- 
act temperature while (fcsT/4J)/ is obtained by fitting the 
spectra to Eq. [S] 



R7v/a(u) cx g v (huj/2)e 



(hQ/2k B T) 



(7) 



In Fig. we show the dependence of Raman spectra 
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with varying characteristic particle numbers at a fixed 
temperature. For small p, most of the weight is located 
in the bulk of the spectrum (u> > 0), whereas, at higher 
p, more and more weight shifts into the tail. For very 
large characteristic particle number, the left edge of the 
spectrum is located at the end of the Co window allowed 
by the support of the density of states and not where the 
Fermi function goes to zero. 




(hco/27i - e°)/4J 



FIG. 5: Raman spectra for four different characteristic parti- 
cle numbers p = 1, 2, 5, 8 at fixed temperature ksT/AJ — 0.1 
in a two-dimensional system. (feT /4J)/ is obtained by fit- 
ting the spectra to Eq. [8] Note that the case where the fit 
does not provide an accurate determination (p = 8) corre- 
sponds to an almost filled band at the trap center, with few 
available thermally excited states. 

From the above descriptions, it is clear that the Ra- 
man spectra strongly depend on temperature and parti- 
cle density. We exploit this strong dependence by fitting 
these Raman spectra to an integrated version of Eq. [5] 

R 2D (Co) t 
* /aK oc - gJM/2) Hi + e (^/2+Mo)/fcflTx (8) 
N p 

From this fit, we extract the temperature and central 
chemical potential of a two-dimensional system. The 
values obtained agree very well with the initial system 
parameters. We summarize the accuracy of the fitted 
parameters in Fig. [6l In the upper panel, we see that, 
asides from small deviations, the temperature can be de- 
termined very accurately by this procedure and is in most 
cases well within 10% of its true value. In the second 
panel, we show that the chemical potential in the center 
of the trap can also be evaluated. Even though its ac- 
curacy is not as good as for temperature, it still agrees 
within 20% for most of the simulated systems. There- 
fore, from the knowledge of experimentally measurable 
Raman spectra, we can accurately determine the tem- 
perature and obtain a good estimate of the chemical po- 
tential at the center of the trap for a free fermionic gas 
confined to a two dimensional optical lattice. Let us point 
out that here we assumed that the hopping amplitude in 
the optical lattice, </, is known. However, as explained 



earlier, J sets the support of the spectra at high fillings 
and can therefore be experimentally detected as a bonus 
by Raman spectroscopy measurements. 
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FIG. 6: Accuracy of detected temperatures and central chem- 
ical potentials. These detected values were obtained by fitting 
each full spectrum with Eq. [8] Upper panel: detected tem- 
peratures are denoted by blue dots while exact temperatures 
by solid black lines. Each shaded region corresponds to a 
10% range centered on the exact temperature and contains 
most fitted points. The regions delimited by dashed lines cor- 
responds to a 20% range. Central panel: central chemical 
potentials are detoned by blue dots, exact values by red a X" . 
Lower panel: density at the center of the trap no. 

In many cases a much simpler fitting procedure can 
already give very good results for the temperature. This 
simplified method focuses on the behavior of the low fre- 
quency part of the Raman spectrum tail. In Fig. the 
Raman spectrum is shown on a logarithmic scale in or- 
der to emphasize its tails. Looking back at the analytical 
expression given by Eq. [5J we see that if e ( huJ / 2 +Ho)/kBT 
is small the Raman signal can be approximated by 

Rl%{Co)*g v {m/2)e^ k ^. (9) 

This simplified expression can be used as long as its valid- 
ity extends over a sufficiently large region of measurable 
signal. In other words, if e ( Su / 2 +c»)/ fc B T <§; l for a wide 
range of HCo/iJ £ [—2; 2]. This condition is most easily 
fulfilled for small or even negative p . Hence, this simpli- 
fied fitting procedure holds best for small and intermedi- 
ate characteristic densities. These limitations are quite 
apparent in Fig. [7] For example, increasing the charac- 
teristic density increases the signal in the tail but drasti- 
cally reduces the region over which the fitting procedure 
works. For p = 5, even at fcsT/4J = 0.1, only a very 
small region is left. The same behavior is also observed 
at larger temperatures. In fact, the good fitting region 
completely drops off the spectrum at large temperatures 
and large characteristic densities as this region would ap- 
pear outside Mj/AJ £ [—2; 2], the range permitted by the 
density of states. The resulting accuracy of the simplified 
fitting procedure is summarized for different system pa- 
rameters in Fig. [5] Here the simplified fit works very well 
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for the lowest values of p whereas for larger characteristic 
densities the temperature is overestimated. 
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p = 8, (k B T/4J) f = 0.560 
p = 5, (k B T/4J) f = 0.121 
p = 2, (k B T/4J) f = 0.096 
p = 1 , (k B T/4J) f = 0.099 
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FIG. 7: Raman spectra in logarithmic scale for four different 
characteristic particle numbers p = 1, 2, 5, 8 at fixed tempera- 
ture fcsT/4J = 0.1 in a two-dimensional system. (fcsT/4J)/ 
is obtained by fitting the tails of the spectra to Eq. [9] 
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FIG. 8: Accuracy of detected temperatures obtained by fitting 
the tail of each spectrum with Eq. [9] Upper panel: detected 
temperatures are denoted by blue dots while exact tempera- 
tures by solid black lines. Each shaded region corresponds to 
a 30% range centered on the exact temperature and contains 
most fitted points. The regions delimited by dashed lines cor- 
responds to a 20% range. Lower panel: density at the center 
of the trap. 



This change in the shape of the Raman spectrum is due 
to the Fermi function whose spread increases with tem- 
perature, fn Fig. [TU] the evolution of the spectra with 
increasing characteristic particle number, p, is shown at 
low temperature. For intermediate filling, the Raman 
spectrum presents its characteristic tail whose size de- 
pends strongly on temperature. At very large filling, the 
tail cannot be followed until its end as the Raman spec- 
trum is limited by the frequency window allowed by the 
cubic density of states. 
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FIG. 9: Raman spectra for three different temperatures 
fesT/6J = 0.1 , 0.5, 1 at a fixed characteristic particle number 
p = 2.5 in three dimensions. (fcsT/6J) 3 is the exact temper- 
ature while (UbT /6J) / is obtained by fitting each spectrum 
up to its peak with Eq. If 01 
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2. Non-interacting fermions in three dimensions 

Let us now look at Raman spectra for atoms confined 
to a cubic lattice. For fixed characteristic particle num- 
ber and various temperatures, typical signals are shown 
in Fig.[9j At low temperatures, these spectra show a main 
peak and a tail that quickly goes to zero. For increasing 
temperatures, the structure of the cubic density of states 
becomes more apparent as weight shifts towards the tail. 



FIG. 10: Raman spectra for three different characteristic par- 
ticle numbers p = 0.5, 2.5, 10 at fixed temperature /cbT/6J = 
0.1 in three dimensions. (fcsT/6J)/ is obtained by fitting 
each spectrum up to its peak with Eq. 1101 

In three dimensions, we can also extract the tempera- 
ture from the measured signal. In this case, the Raman 
rate in the continuum approximation is given by 



g v (K>/2) 
oc 



N 



P 



dp ■ 



(^o-M) 1/2 



.(10) 
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By fitting three dimensional spectra with Eq. [TUl we 
checked that the temperature and central chemical po- 
tential are, to a good degree, uniquely defined for a 
given spectrum. The quality of the extracted temper- 
ature and central chemical potential values are summa- 
rized in Fig. 1111 We find very good agreement between 
the input and extracted temperatures if the fit is done 
from d> m in to the Q value corresponding to the peak of 
the spectrum. However, fitting over the whole spectrum 
is not as successful as more importance is given to the 
rightmost portion of the signal which is not as sensitive 
to temperature as the tail is. Using the reduced fitting 
range, the temperature can be determined within 10% 
uncertainty. The chemical potential at the center of the 
trap can also be determined. However, the agreement 
between the input and extracted values decreases with 
increasing temperature. 
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contrast for larger values of /?, the range over which the 
simplified expression can be fitted becomes very small or 
nonexistent. In Fig. Q31 the extracted temperatures are 
compared to the input temperatures. As expected the 
procedure works well for small and intermediate charac- 
teristic particle numbers. In contrast, for large values of 
p, this method overestimates the system temperature. 




(hco/27i - ep/6J 



FIG. 12: Raman spectra in logarithmic scale for three dif- 
ferent characteristic particle numbers p = 0.5, 2.5, 10 at fixed 
temperature fcsT/6J = 0.1 in three dimensions. (fesT/6J)/ 
is obtained by fitting the tails of the spectra to Eq. [11] 
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FIG. 11: Accuracy of detected temperatures and central 
chemical potentials. These detected values were obtained by 
fitting each spectrum up to its peak with Eq. llQI Upper panel: 
detected temperatures are denoted by blue dots while exact 
temperatures by solid black lines. Each shaded region corre- 
sponds to a 10% range centered on the exact temperature and 
contains most fitted points. The regions delimited by dashed 
lines corresponds to a 20% range. Central panel: detected 
central chemical potentials are detoned by blue dots, exact 
values by red "X" . Lower panel: density at the center of the 
trap. 
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Considering only the tail of the spectrum, a simplified 
fitting procedure can also be used to evaluate the tem- 
perature of a three dimensional gas. Looking back at 
the analytical expression given by Eq. 1 101 we see that 
if e i^/ 2 +^o)/k B T j g sma j] j27j the Raman signal can be 
approximated by 

Rl%{Q)K 9v {t&/2)eW™° T ). (11) 

As in the two dimensional case, we expect this expres- 
sion to be accurate for small or even negative values of 
p a and small Juo/6J 6 [—2; 2]. In Fig. [121 we apply this 
simplified fitting method to spectra with various charac- 
teristic particle numbers. For small values of p — 0.5, 2.5, 
the fit works nicely over a wide range of frequencies. In 



FIG. 13: Accuracy of detected temperatures obtained by fit- 
ting the tail of each spectrum with Eq. 1111 Upper panel: 
detected temperatures are denoted by blue dots while exact 
temperatures by solid black lines. Each shaded region corre- 
sponds to a 30% range centered on the exact temperature and 
contains most fitted points. The regions delimited by dashed 
lines corresponds to a 20% range. Lower panel: density at 
the center of the trap. 

Finally, to conclude this section, we need to point 
out that the frequency resolution attainable in experi- 
ments may not be as good as assumed here. Therefore, 
we checked that the fitting procedure still works for a 
reduced frequency resolution by binning the simulated 
spectra. For two dimensional systems, the temperature 
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extraction method still works surprisingly well. The tem- 
perature can be determined very accurately even if only 
a few points are left on the spectrum (HACj/AJ = 0.8). 
In three dimensions, the fitting procedure still works for 
a frequency resolution of about HAu)/6J = 0.4. 

B. Fermions with (moderate) interactions : 
thermometry from the wings of the density profile 

When atoms in the (|1), |1')) mixture interact via a fi- 
nite interaction strength U, the structure of the Raman 
spectrum changes significantly due to weight redistribu- 
tion in the spectral function (cf. Sec. HVjl. However, up 
to intermediate interaction strengths, the low density re- 
gion on the periphery of the trap is still well described by 
a system of non-interacting fermions. Thus, this region 
can be used to extract the gas temperature assuming the 
system is in thermal equilibrium. The experimental fea- 
sibility of the detection of a small boundary region has 
been shown using radio- frequency spectroscopy for an im- 
balance Fermi mixture and has been used to detect the 
temperature in the absence of an optical lattice potential 
[28| . In contrast to our proposal, the boundary region in 
that case was only occupied by the majority component 
so the Fermi gas in the wings was clearly non-interacting. 
In Fig. Q3] and [TH1 we show density profiles, obtained 
from dynamical mean-field calculations, for two three- 
dimensional interacting systems. We compare these pro- 
files to those calculated by using local density approx- 
imation and the simple Hartree approximation for the 
relation between the density n and chemical potential fi. 
The Hartree approximation simply amounts to inverting 
the relation: [i = ^u=o( rl ) + % n > where ^u=o(^) is the 
chemical potential of the free system for a given density 
n. As one can see from these figures, the interacting 
and Hartree-approximated profiles agree quite well for 
n < 0.3 whereas, for larger densities, the profiles are con- 
siderably different. Hence, by only probing the region 
at the periphery of the trap, we can detect the system 
temperature as the atoms at these locations are still de- 
scribed by a quasi non-interacting model. The validity of 
this approximation will be further evidenced in Sec. IIVI 
Two examples of Raman spectra obtained by collecting 
Raman signal coming from one of the six "semi-spherical" 
regions of low density are shown in Fig. [TBI and [TT1 These 
Raman spectra are simulated using Eq. [5] where the sum 
over positions is limited to one of the six low density re- 
gions and the central chemical potential, /i Q , is the one 
setting the right atom number in the interacting system. 
To show that the temperature can still be well detected 
in this limit, we fit these spectra using the continuum 
Raman expression: 

— ■ cx / / ax sm (bad) 

N P Jo J-$- 

x (121 

A I _|_ e -(hu/2+n a ~x)/(k B T) ■ 



In this expression, where the spatial integral is limited to 
the probed region, there are only three fitting parame- 
ters: the temperature, the central chemical potential and 
a, a parameter related to the size of the probed region 
[291 ] . As shown in Fig. [TBI and [T7l the system temperature 
can be measured successfully using this procedure. 




site index (i,0,0) 

FIG. 14: Density cut through (a;, 0, 0) for an interacting 
system at U/6J = 2, k B T/6J = 0.1 and p = 8.9. Below 
n = 0.3 the interacting and Hartree corrected density profiles 
agree quite well. The shaded region (a; > 30) corresponds to 
the probed area, this region contains 0.4% of the total atom 
number. 




FIG. 15: Raman spectrum obtained by only collecting signal 
from the low density region shown on Fig. [14] The system ex- 
act temperature is ksT/GJ = 0.1 while the detected temper- 
ature obtained through a fit of Ea.[T2"1is {k B T/6J) f = 0.102. 



IV. SPECTRA FOR INTERACTING 
FERMIONS: FROM STRONGLY CORRELATED 
FERMI LIQUIDS TO MOTT INSULATORS 

In this section we discuss the structure of the Ra- 
man spectrum in different strongly-correlated states. We 
present general considerations based on the separate 
spectral contributions of quasiparticle excitations and of 
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site index (i,0,0) 

FIG. 16: Density cut through (a;, 0, 0) for an interacting 
system at U/6J = 2, k B T/QJ = 0.5 and p = 8.9. Below 
n = 0.3 the interacting and Hartree corrected density profiles 
agree quite well. The shaded region [x > 31) corresponds to 
the probed area, this region contains 1.4% of the total atom 
number. 



x10" 3 

(k B T/6J) f =0.517 




FIG. 17: Raman spectrum obtained by only collecting signal 
from the low density region shown on Fig. 1161 The system ex- 
act temperature is k B T/6J = 0.5 while the detected temper- 
ature obtained through a fit of Eq.[PJis {k B T/6J) f = 0.517. 



incoherent high-energy excitations. These considerations 
are illustrated by explicit calculations for the Hubbard 
model with repulsive interactions, treated in the frame- 
work of dynamical mean-field theory (DMFT) The 
DMFT calculations are performed using the numerical 
renormalization group method as an impurity solver [3l| . 
The weakly correlated regime, strongly correlated Fermi 
liquid and Mott insulating regimes are discussed. 

For simplicity, we focus on a homogeneous system, cor- 
responding to Raman spectroscopy being performed in a 
local manner [20( and probing deep inside the bulk of a 
certain quantum state to avoid the influence of a neigh- 
boring state with different character [32j . We also restrict 
our discussion to the paramagnetic phase, and the cal- 
culations are performed at zero temperature (although 
some qualitative remarks will be made on finite temper- 
ature effects). The possible use of Raman spectroscopy 



to detect and investigate the magnetically ordered phase 
is left for future work. 

We focus in this section on zero momentum transfer 
q = 0. This is in contrast to the previous section, in 
which we used q = Tt/a in order to spread the signal as 
much as possible to probe thermally excited states. Here, 
on the contrary, we want to separate and resolve the dif- 
ferent spectral features (e.g. quasiparticles and Hubbard 
bands) as well as possible, and for this q = is more fa- 
vorable [33[ . We focus both on the momentum-resolved 
spectrum (i.e. after time of flight) and on the momentum 
integrated signal (the latter being easier to achieve exper- 
imentally in the lattice) related to the spectral function 
A(k, hu) H| by: 

-R q =o(k, w) = C n F (el - huj + £ k - fi) 

x A(k, el - hw + e k - M)i (13) 

-R q =o( w ) = C J dknp(£3 — hw + £ k — fi) 

x A(k, £3 - Hu + e k - n). (14) 

We recall that, in these expressions, £3 + £k is the dis- 
persion of the outcoupled state |3), while (j, is the chem- 
ical potential of the interacting (|1), |1')) mixture. From 
Eq. Q]the prefactor reads C — 2ir\Q, e \ 2 /h (we note that 
the Wannier matrix element = 1 assuming the 

same lattice potential for state |1) and |3)). 

When specializing to DMFT calculations, the self- 
energy only depends on frequency, so that the spectral 
function 7rA(k, hu) = -Im{l/[^ + ^i-£ k -S(^ + i0 + )]} 
depends on momentum through £ k only. In this case, the 
momentum integration can be replaced by an integra- 
tion over the density of states, g v (s), associated with the 
dispersion £ k (for simplicity, the DMFT calculations pre- 
sented below will be performed for a semi-circular density 
of states): 

i?. q= o(^) = CV J de g v (e)np(el — huj + e — fi) 

x A(k, £3 - hw + £ - fi). (15) 

At T = 0, the Fermi function in these expressions limits 
the integration domain to momenta such that e < fi + 

hu — £3. 

We note that these spectra obey the following sum- 
rules, valid at arbitrary temperature T: 

J dwi? q=0 (k,u;) = Cm(k) , J dujR^uj) = C N/2. 

(16) 

Hence, the total intensity of the momentum-resolved 
signal is proportional to the momentum distribution 
ni(k) = (c k jCic.i) of particles of type |1) in the system, 
while the total intensity of the momentum-integrated sig- 
nal is proportional to the total number of particles N/2 in 
state |1). In these expressions, the frequency integration 
is over the whole range of frequencies where the signal 
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is non-zero (this range is bounded from below, as shown 
later). 

Let us first discuss the shape of the Raman spectrum 
in the simple case of a non-interacting system for which 
A(k, hv) — 8{hv + [i — £k). We obtain in this case: 



C S(hw - £° 3 )n F {e^ - fi), 
N 



(17) 



Hence, at q = and in the absence of interactions, Ra- 
man transitions only exists at the frequency huo = e%. 
This is due to the assumption that the dispersions for 
the atoms in the (|1), |1')) mixture and in the outcou- 
pled state |3) are the same (cf. Fig. [3J . At T = 0, the 
momentum-resolved signal is non-zero only for momenta 
inside the Fermi surface £k < ^ because this spectroscopy 
probes only occupied states. At T ^ 0, the signal extends 
beyond the Fermi surface because of thermal broadening 
according to the Fermi function. These simple consider- 
ations are nicely illustrated by the spectra displayed in 
Figs. US and [19] These results correspond to the Hub- 
bard model with a very low density of particles per site 
(n w 0.18, dilute system) and a very high density of par- 
ticles per site (n w 1.82, or low density of holes in a 
band insulator), respectively. Despite the fact that these 
DMFT calculations were done for a rather high value of 
U/D = 3.5, the system is in effect weakly correlated be- 
cause the density of particles (or holes) is small. Here D 
is the half bandwidth. This is clearly seen from the dis- 
played momentum-resolved spectral functions (Figs. I18[ 
IT51 (b)) as they are weakly modified as compared to the 
non-interacting case. Both spectra show a very sharp 
peak which disperses essentially according to the free dis- 
persion £k (only a shift in position is seen). Hence, the 
Raman spectra are closely following the non-interacting 
behavior: the momentum-integrated spectra ( Figs. [TBI [T9l 
(c)) are sharply peaked, the momentum-resolved Raman 
spectra (Figs. HUES (a)) have very little momentum dis- 
persion (in contrast to the spectral function itself), and 
are suppressed for momenta outside the Fermi surface. 
Due to this narrow momentum dispersion a very sharp 
peak occurs in the momentum-integrated Raman-spectra 
(Figs. HH [19] (c)) (the peak position will be discussed 
later). Let us emphasize that these findings further sup- 
port the detection scheme for interacting particles pre- 
sented in Sec. IIIIBI which relics on the assumption that 
in the low density regions the spectral function behaves 
like the one of non-interacting particles. 

We now turn to spectra in which effects of strong cor- 
relations become more pronounced. In order to discuss 
these spectra on a general basis, we can separate the spec- 
tral function into a contribution from quasiparticles and 
a contribution from high-energy incoherent excitations: 



A(k,hv) =Aj P (k,H + A hlc (k,hv). 



(18) 



The quasiparticle contribution can be appropriately de- 
scribed, at low excitation energies and close to the Fermi 





(a) 














-4 -2 2 4 -4 -2 2 4 



0.6r 



0.8r 



(c) 



(d) 



o 
> 



^? 0.3 



ir 



> 
< 



0.4 



0. L 




o 1 

-4-2 n 2 4 -4 -2 2 4 
(hco/2;i - elp/D hv/(2n D) 

FIG. 18: Spectra for a low density liquid. fi/D — —0.75, 
U/D — 3.5, n ~ 0.18. (a) Momentum-resolved Raman spec- 
trum (in arbitrary units), (b) Momentum-resolved spectral 
function (in arbitrary units), (c) Momentum-integrated Ra- 
man spectrum, (d) Momentum-integrated spectral function. 
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FIG. 19: Spectra for a high density liquid. fi/D — 4.25, 
U/D = 3.5, n ~ 1.82. (a) Momentum-resolved Raman spec- 
trum (in arbitrary units), (b) Momentum-resolved spectral 
function (in arbitrary units), (c) Momentum-integrated Ra- 
man spectrum, (d) Momentum-integrated spectral function. 



surface, by a sharply peaked Lorentzian: 

r k 



A 



qp 



k, hv) ~ — — 

7T [fW 



if 



(19) 



In this expression, e£ p is the dispersion relation of quasi- 
particles, Tk is their inverse lifetime and the spectral 
weight associated with the contribution of quasiparticles 
to the total spectrum of single-particle excitations [35[. 
In a Fermi liquid, the quasiparticle excitations become 
long-lived coherent excitations as the Fermi surface is 
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approached, corresponding to a sharp peak with width 

r k cx( £ £ p - M ) 2 ~(k-k F )\ 

To illustrate how quasiparticles contribute to the Ra- 
man spectrum, we display in Fig. [^D] the results of a 
DMFT calculation for the half-filled Hubbard model at 
U/D = 1.5, which corresponds to a Fermi liquid in the in- 
termediate correlation regime. The momentum-resolved 
spectral function (Fig. [30] (b)) clearly displays a quasi- 
particle peak. This quasiparticle peak becomes sharp 
as the Fermi surface is reached (corresponding here to 
£k F = 0), while for momenta far from the Fermi sur- 
face only a broader incoherent contribution is seen. The 
momentum-resolved Raman spectrum (Fig.[20](a)) shows 
the same features below the Fermi level. However, the 
dispersion of the quasi-particle peak close to kp and 
the incoherent contribution far from the Fermi surface 
behave differently than in the spectral function. The 
momentum-integrated Raman spectrum (Fig. [201(c)) has 
a well-marked peak corresponding to quasiparticle con- 
tribution to the density of states, and a broder hump 
corresponding to incoherent excitations. In order to un- 
derstand better these spectral features, we note that the 
contribution of quasiparticles to the momentum-resolved 
Raman spectrum reads, using (|19p into (|13p: 



^ P =o(k,u;) 



C n F (el P 



r k 



[hiu - 



(20) 



From this expression, it is clear that the quasiparticle 
peak in the momentum-resolved Raman spectrum dis- 
perses according to: (tkj — £3)^ = £k — £ k P ~ (vp — 
Vp ) • (k — kp) + . . . . The last expression is valid for mo- 
menta near the Fermi surface and involves the difference 
between the actual Fermi velocity, v^P, in the presence of 
interactions (related to the effective mass) and the bare 
Fermi velocity, vj?. Indeed, the peak in the Raman signal 
is less dispersive (Fig. |20|) than the one in the spectral 
function (dispersing as vJF • (k — kp)). In practice, since 
the dispersion £ k of the outcoupled band is known, it shall 
be possible to extract directly e^F from the Raman signal 
by plotting it as a function of fku— £3— £ k (as done in [6| in 
the continuum) . As the quasiparticle peak becomes very 
sharp near the Fermi surface k ~ k^, those momenta 
dominate the momentum-integration (at least for lat- 
tices with a non-singular g v (e)). Hence, the momentum- 
integrated spectrum (Fig.[20](c)) has a quasiparticle peak 
located at (fku — s^) qp = (£ kF ) — \i — nu=o — M- I n 
the first expression, (£k F ) corresponds to a Fermi sur- 
face average. The second expression is valid when the 
Fermi surface is only mildly deformed by interactions, so 
that the Luttinger theorem (conservation of Fermi sur- 
face volume) implies that £ kF = ^u=o, where /j.u=o is the 
chemical potential of the non-interacting system at the 
same density. This analysis accounts well for the loca- 
tion of the peak (at ~ —fx) in the spectrum of Fig.l2"01 (c) 
(which corresponds to half-filling, so that [iu=o — 0, 
while fi — U/2 — 0.75 D). The onset of Raman ab- 



sorption in the momentum integrated spectrum at T = 
corresponds to the restriction due to the Fermi function 
hu — £3 > £ k — fi and hence corresponds to a thresh- 
old frequency: (huj — £3) = —D — /i, again well obeyed 
in Fig. [20] (c). We note that this absorption threshold 
corresponds to the transfer of states from the bottom 
of the band with k = (i.e. occupied states well below 
the Fermi surface), as it is known from radio- frequency 
spectroscopy 0,[3- 
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FIG. 20: Spectra for a Fermi liquid with moderate cor- 
relations. fi/D = 0.75, U/D = 1.5, n = 1. (a) 
Momentum-resolved Raman spectrum (in arbitrary units). 

(b) Momentum-resolved spectral function (in arbitrary units) . 

(c) Momentum-integrated Raman spectrum, (d) Momentum- 
integrated spectral function. 



Having discussed a Fermi liquid in the regime of in- 
termediate correlations, we turn to the opposite limit of 
a very strongly correlated system: a Mott insulator, as 
realized e.g. in the Hubbard model at half-filling and for 
large interaction strength (U/D = 3.5 in Fig. l2"Tj) . There, 
quasiparticles are absent and the high-energy incoherent 
excitations correspond to the Hubbard "bands". The 
lower (upper) Hubbard band (LHB, resp. UHB) corre- 
sponds to the process of removing (adding) an atom on 
a singly occupied site. This corresponds to two peaks in 
the spectral function (Fig.l2"TKb) and (d)) at hv^ HB < 
and /u/™ B > 0, separated by the Mott gap A g . Since 
the excitation energy from the ground-state for remov- 
ing a particle is fx, the lower Hubbard band fti/^ HB is 
centered at ~ — fi. This band disperses over a band- 
width of order D with a width r£ HB of order D itself. 
Hence, the excitation is 'incoherent' in nature (except at 
momenta near the top of the band where the width is 
smaller, of order D 2 /U). Similar considerations apply to 
the upper Hubbard band (centered at ~ U — fi). The 
total (momentum-integrated) weight of the lower Hub- 
bard band in the spectral function is proportional to n/2, 
while that of the upper Hubbard band is proportional to 
1 - n/2. 
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At T = 0, the lower Hubbard band is fully visi- 
ble in the Raman spectrum, as seen on Fig. (2D ( a ) 
and (c) (there, -fj, = -U/2 = -1.75 D). This lower 
band in the momentum-resolved Raman spectrum is lo- 
cated at: (fuv - £3)™ = £k - A* - ^k HB - Since 
^lhb < ^ ne i ower Hubbard band is apparent for 
all momenta (in contrast to a quasiparticle peak which 
is suppressed as the Fermi surface is crossed). As the 
momentum integration is dominated by e k — and 
^lhb j g cen t erec [ a t the lower Hubbard band re- 
sults in a peak in the momentum-integrated spectrum 
located at [hw — eg) LHB ~ —fi — (— /x) = 0, as clear 
on Fig. [21] (c). The threshold for Raman absorption 
corresponds to: (tko — e^) th = Min k [e k — ftf LHB ] — fj,. 
The minimum is usually realized for e k = — D, so that: 
(tvjj - 



-o\th 



D 



""top 



fi. At half-filling this reads: 
-D + A g /2 - U/2, with A 9 the Mott 
gap. For temperatures comparable or higher than the 
Mott gap, the upper Hubbard band will become visible 
in Raman spectra, at a location: (hw — e^) 11 ™ ~ — U. 
Concretely, imaging the integrated or full Hubbard bands 
would be very useful as it would not only give informa- 
tion on these incoherent excitations themselves, but as 
well provide a novel method to extract the interacting 
strength, U, and the gap size, A g . 




expected from the analysis above. We note that these 
two features will in general have very different temper- 
ature dependences. As the temperature is raised, the 
quasiparticle peak will be suppressed when temperature 
exceeds the quasiparticle coherence temperature, of or- 
der Zk F D. In contrast, the lower Hubbard band will start 
losing weight (and the upper Hubbard band will start ap- 
pearing) only at a higher temperature scale comparable 
to the gap scale. 

In summary, Raman spectroscopy is a useful probe to 
explore various possible regimes of correlations. Broad 
Hubbard bands are seen in the incompressible Mott 
regime, while the additional observation of a quasipar- 
ticle peak at low temperature signals the formation of a 
strongly correlated Fermi liquid. 
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FIG. 22: Spectra for a strongly correlated Fermi liq- 
uid. n/D = 0.75, U/D = 3.5, n w 0.85. (a) 
Momentum-resolved Raman spectrum (in arbitrary units). 

(b) Momentum-resolved spectral function (in arbitrary units) . 

(c) Momentum-integrated Raman spectrum, (d) Momentum- 
integrated spectral function. 
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FIG. 21: Spectra for a Mott insulator. fi/D = 1.75, U/D = 
3.5, n — 1. (a) Momentum-resolved Raman spectrum (in ar- 
bitrary units), (b) Momentum-resolved spectral function (in 
arbitrary units) . (c) Momentum-integrated Raman spectrum, 
(d) Momentum-integrated spectral function. 

Finally, we display results for a strongly correlated 
Fermi liquid with a spectral function that displays si- 
multaneously a central peak of quasiparticle excitations, 
as well as lower and upper Hubbard bands (Fig. [221 cor- 
responding to a rather large coupling U/D = 3.5 with 
n = 0.85, i.e. to a strongly correlated Fermi liquid). The 
Raman spectra reveal both types of excitations, which 
also lead to two distinct features in the momentum- 
integrated Raman spectrum (Fig. [221 (c)) at frequencies 



V. CONCLUSION 

In this work, we demonstrated that Raman spec- 
troscopy is a versatile probe that can be used to measure 
the temperature of Fermi gases confined to optical lat- 
tices and to identify various signatures of strongly corre- 
lated fermionic phases. The proposed detection scheme, 
implementable with present technology, relies on trans- 
ferring a portion of the atoms stored into the optical lat- 
tice potential to a third hyperfine state. This Raman 
rate can in principle be both resolved in frequency and 
momentum. We showed that momentum resolution is 
not required to accurately measure the temperature of 
free and weakly interacting fermionic atoms loaded into 
an optical lattice and that the detection can either be 
done locally or globally. We also demonstrated that de- 
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tecting several features of strongly correlated liquids and 
Mott insulators such as quasiparticle peaks and Hub- 
bard bands can be done using the same scheme without 
knowledge of the atom momentum. However, in the fu- 
ture, if momentum resolution is experimentally achieved 
in a lattice, momentum resolved Raman rate could pro- 
vide valuable information on the level of correlation of 
fermionic cold atom systems. Finally, we would like to 
point out that Raman spectroscopy can even be used to 
cool down fermionic atoms confined to an optical lattice 
as explained in 36]. 
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